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O ' Generalization of the modified KdV equation to a multi-component system, that is expressed 

^ ■ by ^ + 6 (E^,i=o C,kU,Uk) 1- + ^ = 0, z = 0, 1, • • • , A/ - 1, is studied. We apply a new 

extended version of the inverse scattering method to this system. It is shown that this system 
has an infinite number of conservation laws and multi-soliton solutions. Further, the initial 



value problem of the model is solved. 

> 
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§1. Introduction 

Since the discovery of the inverse scattering method (ISM), it has been shown that the ISM is 
applicable to many soliton equations Among the soliton equations, the modified Korteweg- 
de Vries (mKdV) equation has been studied extensively because of its simplicity and physical 
significance^'BSi^ Generalization of 2 x 2 AKNS formulation to 3 x 3 or, more generally, L x L 
formulation has already been done by a number of authors. i'i'Hi'00'@) 

Generalization of the mKdV equation to a multi-component system or a matrix equation has 
been studied by some authors. One example is a vector version of the mKdV equation 

proposed by Yajima and Oikawa.0^ Sasa and SatsumaS' solved the initial value problem of the 
system, and constructed multi-soliton solution. Another example is a matrix version of the mKdV 
equation studied by Athorne and Fordy.@) 

Recently, Iwao and HirotaEl' discussed a simple coupled version of the modified KdV equation, 

^ + 6 ^^^^^ + ^ = 0, i = 0,l,...,M-l, (1.1) 

\j,k=0 I 

where the constants Cjfc are set to be symmetric with respect to the subscripts, Cjfc = C^j, without 
any loss of generality. They obtained multi-soliton solution of this system with the condition 
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Cjj = 0. HirotallS also studied solutions of semi-discrete version of the model. We call eq. (1.1) 



the coupled modified KdV (cmKdV) equations. The cmKdV equations for M = 1, 2 have been 
solved by the ISM. However, it has not been known whether the cmKdV equations for M > 3 and 
their hierarchy can be solved by the ISM or not. 

In this paper, we propose a matrix generalization of the ISM that includes matrix mKdV equa- 
tion, matrix nonlinear Schrodinger equation (NLS) equation and other integrable equations. This 
formulation also includes the cmKdV equations as a reduction of the matrix mKdV equation. It 
provides us a method to solve the initial value problem and obtain soliton solutions. 

The outline of the paper is as follows. In §2, we introduce a Lax representation for the matrix 
mKdV equation and the matrix NLS equation. Further the reduction of the matrix mKdV equation 
to the cmKdV equations is given. In §3, we perform the ISM for the matrix niKdV equation and 
the matrix NLS equation. In §4, we cast the results in §3 into those for the cmKdV equations. The 
last section, §5, is devoted to the concluding remarks. 

§2. General Formulation 

2. 1 Normalization 

We introduce an M-component vector field u and a constant M x M matrix G, 

u = {uQ,ui,--- ,UM-if , G = {-Cij), (2.1) 

where the symbol T means the transposition. Using this notation, a system of the cmKdV equations 
( pTl] ) is expressed as 

Ut-6{u'^Gu)u^ + u^^^ = 0. (2.2) 

We assume that G is a real symmetric and regular matrix in what follows. Because a real symmetric 
matrix is diagonalized by a real orthogonal matrix, we can put 

G = P^AP, P^P = PP^ = /, (2.3a) 

A = diag(Ao,---,AAf_i), Xj^O. (2.3b) 
Thus, defining a new set of dependent variables v = {vq, vi, ■ ■ ■ , vm-i)^ as 

V = Pu, (2.4) 



Vt-6{v'^Av)v^ + v^^^ = 0, (2.5) 



eq. (p.2D is cast into 



or more explicitly 



1^ + = 0, A,^0, i = 0,l,...,M-l. (2.6) 
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If we change a scale of Vi by \/|A7| • — > Vj, we finally obtain normalized cmKdV equations, 
dt 



i = 0,l,---,M-l. 



(2.7) 



The above equation (2.7) is more convenient than eq. ( |1.1| ) to perform the ISM. Thus, we mainly 
deal with eq. ( p.7| ) as the cmKdV equations in the following. In addition, the dependent variables 
{vi} are assumed to be real. 

2.2 Lax pair 
We consider a set of auxiliary linear equations 



(2.^ 



Here ^ is a {p + (7)-component vector and U, V are {p + q) x {p + q) matrices. The compatibility 
condition of eq. (|2.8D is given by 



Ut-V^ + UV-VU = o. 



(2.9) 



We call U, V Lax pair and eq. (|2.9|) a zero-curvature condition, or Lax equation. We introduce the 
following form of the Lax pair, 



u = K 



-h 










Q ' 




+ 







h _ 




R 






(2.10) 



V = iC^ 



+ 



-4/i O O AQ 

+ c 

O 4/2 iR O 

QxR — QRx —Qxx + 2QRQ 
-Rxx + 2RQR RxQ — RQx 



+ K 



-2QR 2Qx 
-2Rx 2RQ 



(2.11) 



where C is the spectral parameter which does not depend on time, Q = 0. Ii and I2 are respectively 
the p X p and q x q unit matrices ; Q is a p x q matrix (made up of p rows and q columns) ; R is a 
q X p matrix. 

Substituting eqs. (IID and (I2I1I ) into eq. ( |2.9| ), we get a set of matrix equations 



Qt + Qxxx — 3QxRQ — SQRQx — O, 

Rt + Rxxx — 'iRxQR — 2>RQRx = O. 
Suppose that R is connected with the Hermitian conjugate of Q by 

R = eQ\ e = ±l. 



(2.12a) 
(2.12b) 

(2.13) 



3 



Then, eq. ( 2.12| ) is reduced to 

Qt + Qxxx - ^e{Q.Q^Q + QQ^Qx) = O, e 



±1. 



(2.14) 



If we restrict Q to be a real matrix, eq. ( |2.1^ ) becomes equivalent to what Athorne and Fordy 
studied.© We call eq. ( 2.12| ) or eq. ( |2.14| ) the matrix mKdV equation. We consider the ISM for 
eq. ( p.l4| ) with e = -1 in §|. 

2. 3 Lax pair for the matrix NLS equation 
We employ another form of V as 



V = iC^ 



2/i 





+ c 


o 


2Q 








+ i 





2/2 


2R 


O 





-QR Qx 
-Rx RQ 



(2.15) 



Substituting eqs. ( ^loD and ( |2l5D into eq. ( |2.9| ), we get a set of matrix equations 

iQt + Qxx - 2QRQ = O, 
iRt - Rxx + 2RQR = O. 



Under the reduction 



eq. ( |2.16| ) is cast into 



R = eQ\ e = ±l. 



(2.16a) 
(2.16b) 

(2.17) 

(2.18) 



iQt + Qxx-2eQQ^Q = 0, e = ±l. 

We call eq. ( ^.16 ) or eq. ( ^.18| ) the matrix NLS equation.^ By changing the time dependence of 
the ISM in §|3|, we can solve the initial value problem for this system ( ^.18 ) with e = — 1. 



2.4 Conservation laws 

In this subsection, we present a systematic method to construct local conservation laws for the 
matrix mKdV equation and the matrix NLS equation. We start from a special class, p = g = n, of 
eq. m 



^1 

^1 
^2 



U21 U22 
V21 V22 





^1 








" ^1 " 




^2 



(2.19) 
(2.20) 



where all the entries in eqs. (^l9D and ( ^3oD are assumed to be n x n square matrices. The following 
is an extension of the method for the n = 1 case. If we define a square matrix L by 

T = ^2^i\ (2.21) 
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we can prove the following relations from eqs. Q, ( gl9D and ( |23o|) , 

MUuT + ?7ii)}t = {tr(yi2r + Vu)}x, 



U21 + c/22r - r;/: 



rUuT 



(2.22) 
(2.23) 



Equation ( 2.23 ) is interpreted as the matrix Riccati equation. Assuming that U is expressed as eq. 
( p.lOD , we have 

Uu = -iCI, U22 = KI, Ui2 = Q, U2i = R, (2.24) 

where Q and R are square matrices in this case. Then eqs. (|]2|) and ( ^33|) are cast into the 
following equations, 

{tr {QT)}t = {some function of T, Q, R and C}x, (2.25) 

2iCgr = -QR + Q{Q-^ ■ QT)^ + (Qr)^. (2.26) 

Equation ( |2.25| ) is nothing but a local conservation law. It shows that tT{QT) is a generating 
function of the conserved densities. We expand QT with respect to the spectral parameter Q as 
follows, 



Substituting eq. ( 2.27| ) for eq. ( p.26| ), we obtain a recursion formula. 



i-i 



Fi+i = -6i,oQR + Q{Q-^Fi), + FkFi^k, 1 = 0,1,- 



(2.27) 



(2.28) 



fc=i 



Each ti Fi is a conserved density for all positive integer /. Using the above formula ( |2.28| ), first four 
conserved densities are given by 

Fi = -QR, (2.29) 

trFa = tr{-QR^}, (2.30) 

trFg = tT{-QR^^ + QRQR}, (2.31) 

trF4 = tr{-QR^^^ + 2QR^QR + QRQ^R + 2QRQR^}. (2.32) 

It should be noted that all elements of Fi = —QR are conserved densities for the matrix mKdV 
equation ( p. 12 ) and the matrix NLS equation (2.16). This fact can be proved simply by a direct 
calculation. 

2.5 Hamiltonian structure of the matrix mKdV equation and the matrix NLS equation 
In this subsection, we consider the matrix mKdV equation and the matrix NLS equation with 



the condition that Q and R are n x n square matrices. The matrix mKdV equation (2.12) has the 
following Hamiltonian structure. A set of the Hamiltonian and the Poisson bracket is 



H = tiJ {iF4}dx = trj l^-iQRccxx + i^QR{QRx - Qx-R)|dx, 



(2.33) 
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and 



where {X Y}f 



{Q{x)fQ{y)} = {R{x)fR{y)} = 0, 
{Q{x)(fR{y)} = i6{x-y)U, 
{Xij,Yki} for matrices X, Y and IT denotes an 



!ki — i^'-tj 

For the matrix NLS equation, the Hamiltonian is 



H' = tv J {-F2}dx = ti J {QRxx - QRQRjda 



(2.34a) 
(2.34b) 

permutation matrix. 



(2.35) 



instead of eq. ( |2.33 ), while the Poisson bracket is the same. We can rewrite Poisson bracket between 
each element of Q and R explicitly as follows, 

{Qij{x),Qki{y)} = {R^j{x),RM{y)} = 0, 



{Qij{x), Rki{y)} = i6ii6jkS{x - y) = in^^/5(x - y). 



(2.36a) 
(2.36b) 

Indeed, we can check the following relations, 

Qt = {Q, H} = -Q^^^ + SQ^RQ + SQRQx, (2.37) 

Rt = {R, H} = —Rxxx + 3RxQR + SRQRx, (2.38) 
which are nothing but eq. ( |2.12| ). It shows that eq. ( 2.33| ) and eq. (2.34) are respectively the 



Hamiltonian and the Poisson bracket for the matrix mKdV equation. The same thing is true for 
the matrix NLS equation. 

2.6 r-matrix representation of the matrix mKdV equation and the matrix NLS equation 

In §2.4, we have shown that the matrix mKdV equation ( 2.12| ) and the matrix NLS equation 
( p. 16 ) have an infinite number of conservation laws. In this subsection, we show that all the integrals 
of motion are in involution. In the following, we consider the systems with an infinite interval and 
assume the rapidly decreasing boundary conditions. 



Q{x,t), R{x,t) — > O as X — > ±oo. 
If we define a classical r-matrix by 

n 

r(Ci,C2) 



(2.39) 



2(Ci - C2) 



o n 
n o 



n 



(2.40) 



the following relation 
{U{x-Ci)^U{y;C2)} = S{x-y) 



r(Ci,C2), U{x;Ci)i 



I 







/ 





®U{x;C2) 






+ 









I 







I 



(2.41) 
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is satisfied. Here [A, B] = AB — BA is a commutator and U is given by eq. ( |2.10| ). By means of 
eq. ( |2.41| ), we obtain a relation for transition matrices,© 



{T{x, y; Ci) f Tix, y; C2)} = [r(Ci, C2), T(x, y; Ci) T{x, y; (2)] 
where T(x, y; C) is the transition matrix defined by 



r(x,y;C) = Pexp|^V(z,C)dz|, 



(2.42) 



(2.43) 



with P being the path ordering. 



Taking the traces of both sides of eq. ( 2.42 ) , we get 

{logr(Ci),logr(C2)} = 0, 

where r(C) is defined by 



(2.44) 



t(C) =trr(oo,-^;C), (2.45) 
for the systems with an infinite interval. Expanding ( 2.44| ) with respect to the spectral parameters 



Ci and C2, we have the involutiveness of conserved quantities {Jn}, 



(2.46) 



This fact indicates the complete integrability of the matrix mKdV equation and the matrix NLS 
equation. 

2.7 Reduction of the Lax pair for the coupled modified KdV equations 

In this subsection, we show a method to reduce the matrix mKdV equation to the cmKdV 
equations. We recursively define 2™~^ x 2™~^ matrices Q^™) and i?^'") by 



q{i) = ^qvq + ivi, R^^^ = eii^iQVo - ivi), 

Q{m) 



(2.47) 
(2.48) 



(2.49) 



-{fJ'2mV2m — i^^2m+l )-^2™--l 

—£2m+l{fJ'2mV2m + ™2m+l)l2"^-'^ 
-{li-2mV2m " if 2m+l )-f2™-l -Q^™) 

Here J2m-i is the 2™~^ x 2™""^ unit matrix. Each ei is a constant which is either 1 or —1. //2m is 
a constant that satisfies 

(2.50) 



2 ^2im 

A*2m ~ T ~ e2mS2m+l- 



£2m+l 



For g(™) and E^*") defined by eqs. (^)- (|2lg|) , we can prove a simple relation. 



2m- 1 



Q{m)^{m) ^ rMqM ^ ^.^2 . 

j=0 



(2.51) 
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by the induction method. Then substituting Q^"^^ and R^"^^ for Q and R in the matrix mKdV 
equation (|2.12|), we obtain 




±1, 



0,1, 



,M-h 



(2.52) 



where we set M = 2m. 
In the following, we choose 



Ei = —1, i = 0,1, ■ ■ ■ , 2m — 1, 
/^2i = 1, j = 0, l,---,m- 1, 



and consider a self-focusing type of the cmKdV equations (|2.7[) , that is, 



dt 



+ 6 E 



+ 



,i=o 



dx dx^ 



0, 



0,1,---,M- 1. 



In this case, the recursion relations of Q^"^) and R^'^^ are 

Qi^)=VQ + wi, R^^^ = -VQ + ivi 



Q(m+1) 
^(m+1) 



1^2™ + i^^2m+l)-^2'"-i 



Vo + \Vi, 
Q{ra) 

{V2rn - i-y2m+l)-^2™-i -Q^"") 

We can show that a simple relation between Q^"^^ and i?^'") holds. 



(2.53) 
(2.54) 

(2.55) 

(2.56) 
(2.57) 

(2.58) 
(2.59) 



It should be noted that the Hamiltonian and the Poisson bracket for the matrix mKdV equation 
( ^.12 ) become invalid for the cmKdV equations ( ^.55| ) with M > 3. This is because the degree 
of freedom of Q^*") and R^'^^ for the cmKdV equations is less than that for the matrix mKdV 
equation. 

§3. Inverse Scattering Method 

In this section we consider the scattering problem associated with 2n x 2n matrix ( 2.10| ) under 
the constraint ( p. 13 ) with e = —1 and the boundary conditions (2.39), that is, 

" -Ki Q 

R iCI 

Q,R{=-Q^)-^0 as x^±oo. (3.2) 

The results of this section are applicable to the matrix mKdV equation, the matrix NLS equation 
and other members of the hierarchy only if Q and R are n x n square matrices. The main idea in 
what follows is a modification of the analysis in refs. 2C for the matrix KdV equation. 



r = -q\ 

±oo. 



(3.1) 
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3.1 Scattering problem 



Let ^'(C) and <&(C) be solutions of eq. (3J) composed by 2n rows and n columns. We can show 
that 



dx 



{^!\C)H^)} = o. 



(3.3) 



Hence we employ the following definition of x-independent matrix Wronskian <1>], 

W[^,^ = ^\C)^{C). (3.4) 
We introduce Jost functions 0, and ^, tp which satisfy the boundary conditions, 



and 



I 
o 

o 



o 
I 

I 
o 



e as X 



as X 



-oo, 



-oo, 



e ^ as X 



+ 00, 



-\C,x 



as X ^ +00. 



(3.5a) 
(3.5b) 

(3.5c) 
(3.5d) 



Here O and / are respectively the n x n zero matrix and the n x n unit matrix. It can be shown 
that (j)e^'^^, 'ipe~^^^ are analytic in the upper half plane of (, and 0e~'^^, tpe^^^ are analytic in the 
lower half plane of ( because Q, R are assumed to go to O sufficiently rapidly at x — > ±00. We 
assume the following integral representation of the Jost functions %[) and 



o 
I 



/•CO 

e'^^+ / K{x,s)e'^'ds, 



I 

o 



e-iC^+ / i?(x,s)e-'^Ms, 



(3.6) 
(3.7) 



where K{x, s) and K{x, s) are column vectors whose elements are n x n square matrices, 



K{x,s) 



Ki{x,s) 
K2{x,s) 



,K{x,s) 



Ki{x,s) 
K2{x,s) 



(3. 



We substitute eq. ( p^ ) for eq. ( p7l| ) and get the relations for Ki and K2, 



lim 


Ki{x,s) 











s— >+oo 


K2ix,s) 








2Ki{x,x) = Q{x), 



(3.9) 
(3.10) 
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{dn, - ds)Ki{x,s) = Q{x)K2{x,s) {s > x), 
{d^ + ds)K2{x, s) = R{x)Ki{x, s) {s > x). 
Similarly, substituting eq. ( |3.7D for eq. (|3.lD, we get for Ki and K2, 

lim 



Ki{x,s) 







K2{x,s) 








-2K2{x,x) =R{x), 
(d, - ds)K2{x, s) = R{x)Ki{x, s) {s>x), 
{d^ + ds)Ki{x,s) = Q{x)K2{x,s) {s> x). 



(3.11) 
(3.12) 

(3.13) 

(3.14) 
(3.15) 
(3.16) 



Because a pair of the Jost functions (p and (p, or ip and ip forms a fundamental system of solutions 
of eq. (|3.l|) , we can set 

cPix, C) = ij{x, C)A{C) + ^Pix, C)B{0, (3.17a) 

^{x, C) = ij{x, OB{C) - ^{x, C)i(C). (3.17b) 

Here the coefficients ^(C)) ^(0) -^(C) -B(C) are x-independent n x n matrices and called 
scattering data. 

According to the asymptotic behaviors of the Jost functions ( |3.5a| )-( |3.5d ), we get 

W[<p, (/>] = W[^, 0] = W[ip, ip] = W[^p, ip] = I, (3.18a) 

W[^, (p] = W['^, V'] = O, (3.18b) 

A{C) = W[i;,<P], (3.18c) 

A{0 = -W[i;,4>], (3.18d) 

B{C) = W[i;,4>], (3.18e) 

BiC) = W[ij,^]. (3.18f) 

The expressions ( 3.18cD -( p.l8d ) show that A{() and ^(C) are analytic respectively in the upper 
half plane and in the lower half plane. Using the above relations ( ^.18a )-( 3.18f ), we obtain the 
following relations among A{(), A(C), B{() and B{C), 



A^C)AiO + B^C)BiC) = I, 
A^{aA{C) + B\C)B{C) = I, 

a\C)b{c)-bHC)a{C) = o. 



These relations are written as 

aHC) 
bHC) 



bHC) 
-aHC) 



Aio m 

B{C) -AiO 





I 










I 



(3.19a) 

(3.19b) 
(3.19c) 

(3.20) 
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which leads to the inversion of eq. ( p.l7| ), 

i,{x, C) = OA\C) + C)^^(C*), (3.21a) 
V'Cx, C) = <P{x, C)bHC) - Hx, OA^C)- (3.21b) 



3.2 Gel'fand-Levitan-Marchenko equation 

To derive the formula of the ISM concisely, we assume that j4(C), A{Q, B{Q and -B(C) are entire 
functions. This assumption is true if the potentials Q and R decrease faster than any exponential 
function at x — > ±00. The result is, however, applicable to larger classes of potentials Q and R. 

Multiplying A(()~^ and A{()^^ from the right to eqs. ( 3.17a| ) and ( p. 17b ) respectively, we get 



</.(x,C)A(C)-' 

4>ixx)Aicr' 



i;{x,C)+i^{xX)B{C)A{0-\ 

-V(x,C) + V^(x,C)^(C)i(C)"'. 



We operate 



1 

2^ 



dCe 



(y > x) 



(3.22a) 
(3.22b) 

(3.23) 



to eq. (|3.22aD , where C is a semi-circle contour from — 00 + 10"*" to +00 + 10"'' passing above all poles 
of {detA(C)}~^. After a standard calculation, we get the Gel'fand-Levitan-Marchenko equation. 



K{x,y) + 
where F{x) is defined by 



O 
I 



F{x + y)+ / K{x,s)F{s + y)ds 



O 
O 



{y > x), 



F{x) = ^j^dCe<-B{C)A{0-\ 
We remark that A{C,)~^ is given by 

1 



:^(C), 



(3.24) 



(3.25) 



(3.26) 



det^(C)' 

where A is the cofactor matrix of A. We assume that l/detj4(C) has N isolated simple poles 
{Ci; C2; • • • ) Ca^} in the upper half plane and is regular on the real axis. Each of these poles determines 
one bound state. Then, by the use of the residue theorem, we get an alternative expression of F, 

N 

(3.27) 



F{x) = — / dCe'^-B{C)A{0-' -iY.Cje"^'''- 



Here Cj is the residue matrix of B{C,)A{C,) ^ at = Cj- 
Similarly, we operate 

" dCe"''^*' (y > x) 



1 

2^ 



c 



(3.28) 
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to eq. (|02b| ), where C is a semi-circle contour from —00 + iO to+oo + iO passing below all poles 
of {detA(C)}^^. We get the counterpart of the Gel'fand-Levitan-Marchenko equation, 



where F{x) is defined by 



O 



_ POO _ _ 

F{x + y)- K{x,s)F{s + y)ds 

Jx 



O 
O 



{y > 3:), 



F{x) = — j_dCe--<-B{Om-'- 



(3.29) 



(3.30) 



If we assume that l/detA(C) has isolated simple poles {Ci, C2; • ■ ■ , Cn} the lower half plane 
and is regular on the real axis, we get an alternative expression of F, 



_ N _ 

Here is the residue matrix of B{QA{(^)~^ at = ^f^. 



(3.31) 



3.3 Time- dependence of the scattering data 

Under the rapidly decreasing boundary conditions ( |3.2D , the asymptotic form of the Lax matrix 
V is given by 

-AiC^I O 
O 41(3/ 
We define time-dependent Jost functions by 



V 



as X ±00. 



(3.32) 



= 0e 



= 0e 



-4iC^t 



o 
o 



o 
I 

I 
o 



as X — > —00, 



e^«^+^^^'* as x^+00, 



as X ^ +00. 



From the relations 



we get 



dt 



dt 



We substitute the definitions of the scattering data, 

(t>{x, C) = i^ix, C)A{C, t) + ^(x, C)B{C, t), 



(3.33a) 
(3.33b) 
(3.33c) 
(3.33d) 

(3.34) 
(3.35) 

(3.36a) 
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4>{x, C) = C)B{C t) - ^{x, C)A{C, t), (3.36b) 
into eq. ( p. 35 ). Then taking the hmit x +00, we get 

At{C,t) = 0, 

Bt{C,t) = 8ieBiC,t), 

and 

MC,t) = o, 

Bt{C,t) = -8iC'B{C,t). 
The above relations lead to the following time-dependence of the scattering data: 

A{C,t)=A{C,0), 

B{C,t) = B{C,0)e^-<'\ 

and 

i(C,t)=i(C,o), 

B{C,t) = BiC,0)e-''<''. 

By eqs. ( |3.38a )-( 3.38d ), the time-dependences of BA^^, Cj and BA^^, Ck are respectively given 

by 



(3.37a) 
(3.37b) 

(3.37c) 
(3.37d) 

(3.38a) 
(3.38b) 

(3.38c) 
(3.38d) 



B{i,t)A{^,t)-^ = B{i,^)A{i,Q)-^e^'^'\ 
C,it) = CjiO)e'"^'\ 



and 



m,t)A{^,t)-' = S(e,0)i(^,0)-ie-«'«'*. 

To summarize, we obtain explicit time-dependent forms of F{x,t) and F{x,t), 

F{x,t) = — r deei«^+8i«'*S(e,0)^(e,0)-i -iV C,•(0)e'^^^■+^^'^^ 
2tt J -00 



F{x,t) = — d^e 
^vr J -00 



N 



5(e,0)^(^,0)-i+i5](7fc(0)e-'' 



(3.39) 
(3.40) 

(3.41) 
(3.42) 

(3.43) 
(3.44) 



k=l 
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3-4 Initial value problem 

Because of the constraint R = —Q^, we have some relations which make the further analysis 
simple. 

First, we have 

det^(C) = {detA(C*)r, (3.45) 

which is proved in Appendix A. This relation gives us a useful information about the total number 
and the positions of the poles of A(()^^ and A{()^^, 

N = N, Ck = Ck- (3-46) 

Second, due to eq. ( p. 19c ), we have 

mMcr' = {B{C)Ain-'}\ (3.47) 

which leads to 



BiOMCr' = {B{OA{Cr'}^ (e : real), (3.48a) 
Ck = Cl. (3.48b) 
The relations j^A^ and (|1|) give a connection between F{x,t) and F{x,t), 

F{x,t) = F{x,t)K (3.49) 

Combining the above results, we arrive at 

poo roo 

Ki{x,y;t) = F{x + y,t)^ - dsi / ds2Ki{x, S2;t)F{s2 + si,t)F{si + y,t)^ , (3.50) 

Jx J X 

K2{x,y;t) = -F{x + y,t) - dsi / ds2K2ix, S2]t)F{s2 + si,t)^ F{si + y,t), (3.51) 

where F{x,t) is given by eq. ( |3.43| ). 

We can solve the initial value problem of the matrix mKdV equation as follows. 
(1) For given potentials at t = 0, Q{x,0) and R{x,0) which satisfy R{x,0) = —Q{x,0)^, we solve 
the scattering problem (3.1), and obtain scattering data {B{^)A{^)~^ XjACj}■ 
{2) The time-dependence of the scattering data is determined by eqs. ( |3.39| ) and ( p. 401 ). 
(3) We substitute the time-dependent scattering data into the Gel'fand-Levitan-Marchenko equa- 
tions ( 3.5C| ) and ( 3.51 ). Solving the equations, we reconstruct the time-dependent potentials, 

Q{x,t) = -2Ki{x,x;t), (3.52) 

R{x,t) = -2K2{x,x;t). (3.53) 

In this way, we obtain the solution Q{x, t) and R{x, t) from the initial condition Q{x, 0) and R{x, 0). 
This procedure proves directly the complete integrability of the matrix mKdV equation (|2.14| ) with 
e = -1. 
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If we employ other time-dependences of the scattering data, for instance, 

BiC, t)A{t tr' = BiC, 0)AiC, 0)-'e^'i'\ (3.54) 

Cj{t) = Cj{0)e'^'^^\ (3.55) 

the initial value problem of the matrix NLS equation ( |2.18 ) with e = —1 can be solved. 

As for the constraint R = —Q^, one comment is in order. Because F is connected with F by eq. 
( p. 49 ), we can prove by the Neumann-Liouville expansion (see Appendix B) that the solution of 
eqs. ( 3.50|) and (|3.51|) satisfies 

K2{x,x;t) = —Ki{x,x] t)^. (3.56) 
This relation assures that the relation R{x,t) = —Q{x,t)'^ holds at any time t. 

3.5 Soliton solutions 
Assuming the reflection-free potentials which satisfy 

B{i) = B{0 = (e : real), (3.57) 

we can construct soliton solutions of the matrix mKdV equation. In this case F{x) is given by 

(3.58) 



N 



F{x,t) = -\Y,C,{t)e<^\ 
i=i 

Cj{t) = Cj(0)e^"^^. 
To solve eq. ( p. 501 ) with eq. ( |3.5^ ), we set 



N 



(3.59) 



(3.60) 



k=l 



Introducing eq. ( |3.60| ) into eq. ( |3.5[lD , we have a set of algebraic equations, 



N N 



p,ix,t)-j2J2 



1=1 j=i (0 Cfc)(Ci cr) 



We define a matrix S by 



N 



Sik = Siki - 



,ti(o-a)(o-cr) 



Pz(x,t)Q(t)tC,(t)e2'(0-C)- = /. 



-Ci{t)^Cj{t), l<l,k<N. 



(3.61) 



(3.62) 



Then eq. ( 3.61| ) is expressed as 



{P1P2 ■■■ Pn 



( S 



11 



\ Sni 



Sin ^ 
Snn ) 



(3.63) 
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Thus the A^-soHton solution of the matrix mKdV equation (p.l4|) with e = — 1 is given by 



Q{x,t) = -2Ki{x,x;t) 

N 

= -2i5]Pfc(x,t)C,(t)te-2iC: 



k=l 



I Ci(t)te-2icr^' ^^ 

C2(t)"^e-2iC2*^ 



(3.64) 



As a special case = 1, we get 1-soliton solution of the matrix mKdV equation, 

„8i{Ci^-Ci*^)t 1 ~^ 

Q{x,t) = -2\l I - — -— Ci(0)tC7i(0)e2'(^i-^i*)" ) Ci{<d)U-^'^*^''-^"^'^ ' 



(Ci-cr)2 

-2i|e-^('^i^'^i*)''-^^(^?-'^i')*/ 
■C7i(0)te-'(^i+^i)^-^'(^?+^i*')*. 



} 



j^^^ Ci (0) t Ci (o)ei(^i -fr)-+4i(c? -cr 



(3.65) 



If we replace the time dependence ( |3.59| ) in eq. ( 3.5S| ) with eq. ( 3.55| ), we obtain the A^-soliton 
solution of the matrix NLS equation ( |2.18| ) with e = —1. 

§4. Reduction of the ISM for the Coupled Modified KdV Equations 

In order to make the ISM in §^ applicable to the cmKdV equations, we have to take into account 
the internal symmetry of Q and R defined by eqs. (|2.56|) -( p.58|) . 
If we set Q^*") and for m > 2 as 

2m-l 2m-l 

g(-) =vol+ ^kek, R^""^ = -vol+ E ^kek, (4.1) 



k=l 



k=l 



the following important relations for 2™ x 2™ matrices {ej} hold 



tr efc = 0. 



(4.2) 
(4.3) 
(4.4) 



Here {■,•}+ denotes the anti-commutator. I is the 2™ ^ x 2™ ^ unit matrix. Equation ( [4.2| ) leads 
to 

tr(e,e,) = -2^-^6ij. (4.5) 



The results in assures that the cmKdV equations have an infinite number of conservation 
laws. From the explicit forms of the conserved quantities, we find that the first four conserved 
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densities for the original cmKdV equations (|1.1| ) are given by 

^CjkUjUk, (4.6) 

UjUk,x, Vj, k (j / k), (4.7) 

^CjkUjUk\ -^CjkUj,xUk,x, (4.8) 

) j,k 

CjkUjUk - 3 ^ CjkUjUk ■ ^ CjkUj^xUk,x + r) ^ ^ Cjk'Uj^xx'^k,XX r) 
\3,k J j,k j,k j,k / x) 

(4.9) 

We remark that the method in §2^ does not give the quantity ( [4.71 ). 

Next, we discuss the initial value problem and the soliton solutions of the cmKdV equations. 
Considering the scattering problem {3A_) with the potentials Q^"^^ and R^"^^ for m > 2, we can 
show that there are following restrictions on the scattering data. 
Proposition 4.1 

(1) The determinant of ^(C) satisfies 

detA(C) = {detA(-C*)}*, (4.10) 

as a function of complex Thus the poles of l/detA(C) in the upper half plane should appear on 
the imaginary axis or as pairs which are situated symmetric with respect to the imaginary axis. 
Therefore, we can set the values of 2N poles as 

C2j~i = Cj+^Vj, j = l,2,---,N, (4.11a) 
C2j = -C2J-1 = -^1 +^Vj, j = 1, 2, • • • , iV, (4.11b) 

for r]j > 0. The condition ( [4.11] ) should be interpreted as follows; if Ci is pure imaginary, it does 
not need its counterpart. 

(2) The reflection coefficient B{^)A{^)^^ for real ^ should be expressed as 

2m-l 

B(^)^(^)-i =rWl+ J2 r^^'^^k. (4.12) 

k=l 

Here r^^^ and are complex functions of ^, t which satisfy 

rW(-^) = rW(0*, r(^)(-0 = rW(0*. (4.13) 

(3) The residue matrices {Ci, C2, • • • , C2N-1, C2n} should be expressed as 

2m- 1 

iC2,-i = cf I + J2 J = 1, 2, • • • , AT, (4.14a) 

k=l 
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2m- 1 

iC2, = cf *I + ^ cf ) *efe, j = l,2,...,N, 



(4.14b) 



fe=i 



where c^°\ cf^ are complex constants. For example, {Ci, C2, ■ ■ ■ , C2iv} for the 4-component 



cmKdV equations are given by 



iC2j-i 



2j 



(4.15a) 
(4.15b) 



in a different notation. A proof of the statements is given in Appendix C. 
Considering the above conditions, we have explicit expressions of F and F in terms of I and e^, 

2JV 

(4.16) 



F{x,t) = — / B{i,t)A{£,,t)-^e^^-di-iY,Cj{t)e<^- 

^TT J -00 

1 roo ( 2m-l 

= _ y j (^(0)gi«:c ^ ^(0) J ^ (^(fc)ei€^ + ^(fc) *e-i«^)efc } d^ 



-ES (4"^e'^^- + cf V'^.^)I+ ^ (cf e'^^- + 

i=l I A;=l 



c}"' e -.-)efe ^, (4.17) 



F(x,t) = F{x,t)^ 

= - / 5(e,i)^(C,i)-'e-«Me + iEC'.(*)^e-<.-' 



2iV 



(o)^iS^ ^ ^(0) *e-i«a;) J _ ^ (^l^ig'^^ + ^l^i *e~'^'^)efc )> d^ 

k=l 



2m- 1 



(r^^^e 



1 

2^ Jo 

Af f 2m- 1 

-E (cf e^^^- + cf V^^^)I- ^ (cf e^O- + cf *e-^^)e, ^. 



Because B{^,t)A{^,t) ^ and Cj(i) depend on t as 



the time-dependences of r^^^ r^'^^ and Cj'\ c"p are given by 

r(0)(^,t) = r-W(^,0)e8'«'*, rW(e,t) = r('=)(e, 0)e8if' 
cf(t) = cf(0)e«'^l*, cf)(t)=cf)(0)e«^^. 



(4.18) 



(4.19) 



(4.20) 
(4.21) 

(4.22) 
(4.23) 
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It should be noted that t) and F{x, t) are expressed as 



2m- 1 



k=l 
2m- 1 



fc=i 

where the real functions f^^\x,t) and satisfy 



(4.24) 
(4.25) 

(4.26) 



Taking into account the conditions ( [4.1l| )-( [l.l4 ), we can advance the analysis in parallel with the 
discussion in §^ for the matrix niKdV equation. Thus the initial value problem of the cmKdV 
equations can be solved by the ISM, as has been shown in §^. 

We replace in §^ by 2N and find the A^-soliton solution of the cmKdV equations (2.55) with 
M = 2m components, 

Q(™)(x,t) = -2Ki{x,x;t) 

2N 

= -2i5]Pfc(x,t)Cfe(t)te-2i«^ 



k=l 



-2i{I I ■■■ I) S-^ 

2N 



I Ci(t)te-2if^ \ 
C2(t)^e-2iC2*^ 

^ C2^(t)te-2iaiV- j 



(4.27) 



where the matrix S is given by 



2N 



Sik = ^iki — 



1 (0-0(0-0 



Ci{t)^Cj{t), l<l,k<2N. 



(4.28) 



It is not evident whether eq. ( 4.27| ) can be expressed as 

2m-l 

Q^"'\x,t) = voix,t)l+ ^ Vk{x,t)ek, 



(4.29) 



k=l 



without using SiCj, eiejCk, etc. But, noting the fact that summations and products of real quater- 
nions are real quaternions, this can be proved for m = 2 (4-component cmKdV equations) by 
using Neumann-Liouville expansion (see Appendix B). It is an open problem to prove eq. (|4.29 ) 
for general M = 2m. 

§5. Concluding Remarks 

In this paper, we have constructed an extension of the ISM to solve the matrix mKdV equation 
and the matrix NLS equation. We get the coupled mKdV equations ( p.7| ) as a reduction of the 
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matrix mKdV equation. Through the extension, we have shown that the coupled niKdV equations 
have an infinite number of conservation laws and multi-soliton solutions and that its initial value 
problem is solvable. 

The existence of conserved quantities for the coupled mKdV equations has been proved by Svi- 
nolupov in a different approach.^) Iwao and Hirota obtained Pfaffian representation of A^-soliton 
solution for the model by means of so-called Hirota's method.© We stress that the initial value 
problem of the coupled mKdV equations has been solved in the present paper for the first time. 
In addition, it directly proves the complete integrability of the model. Our scheme enables us to 
construct more general solution than the already known solutions. 

We can transform the coupled mKdV equations to a new integrable coupled version of the Hirota 
equation, ^ by changes of dependent variables vj and independent variables x, t. The new coupled 
Hirota equations describe interactions among different modes in optical fibers and seem to be 
physically significant. Wide applicability of our extension of the ISM will be reported in the 
succeeding papers. Further, the problem of integrable boundaries at x = for the coupled mKdV 
equationJll^ will be studied elsewhere. 

After completing writing the paper, the authors were informed by Hisakado that the similar Lax 
formulation was used in the work of Eichenherr and Pohlmeyer.§l) 
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Appendix A: Proof of eq. ( 3.45 ) 



We begin with a counterpart of the linear equations (|2.8|) , 

= U^, (A-1) 

where ^' is assumed to be a square matrix. Using this equation, we get a chain of identities, 

*^^'-^ = ;7, tr^og^-)^ = trC/, 
(logdet^-)^ = trC/, 

det^- = det*(xo) • e*''^-o ^'^'^ . (A-2) 
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In case that U is given by 

-iCI Q 
R iCI 

with square matrices /, Q and R, eq. (|A-2|) leads to 



U 



det^ = const. 



If we take [ </> ^/^ ] as ^, we get 



that means, 



det[ V ] = det^^(C*) = det^(C), 



(A-3) 
(A-4) 

(A-5) 
(A-6) 



detA(C) = {det^(C*)}*. 
Appendix B: Neumann-Liouville Expansion 

Due to the Gel'fand-Levitan-Marchenko equations (|3.24|) and (|3.29D , we get closed integral equa- 
tions for Ki, K2 and 7^2, Ki, 

_ roo roo _ 

Ki{x,y)=F{x + y)- dsi ds2Ki{x, S2)F{s2 + si)F{si + y), (B-1) 

Jx Jx 

_ roo roc _ _ 

K2{x,y) = -F{x + y)- dsi ds2K2{x, S2)F{s2 + si)F{si + y), (B-2) 

Jx J X 

roo _ roo roo _ 

K2{x,y) = - dsF{x + s)F{s + y)- dsi ds2K2{x, S2)F{s2 + si)F{si + y), (B-3) 

Jx J X J X 

_ roo _ roo roo _ _ 

Ki{x,y) = - dsF{x + s)F{s + y)- dsi ds2Ki{x, S2)F{s2 + si)F{si + y). (B-4) 



By successive approximations, we obtain the Neumann-Liouville expansions for Ki, K2 and K2, 



Ki{x,x) = F{2x) 



dsi ds2F{x + S2)F{s2 + si)F{si+x) 



+ ■■■ 

H-'i-T I dsi / ds2- 



dS2nF{x + S2n)F{s2n + S2n-l)F {S2n-l + S2n-2) 

■■F{s2 + si)F{si + x) 



+ ■ 



(B-5) 



K2{x,x) = -F{2x) 



dsi ds2F{x + S2)F{s2 + si)F{si+x) 
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roo poo roc _ 

+ / dSi / dS2---/ dS2nF{x + S2n)F{s2n + S2n-l)F{s2n-l + S2n-2) 

Jx J X J X 

■■■F{s2 + si)F{si+x) 
+ ■■■, (B-6) 

roo _ 

K2{x, x) = - dsiF{x + si)F(si + x) 

J X 

poo poo poo _ _ 

+ dsi dS2 dS3F{x + S3)F{S3 + S2)F{S2 + si)F{si + x) 

J X J X Jx 



IX J X 

+ ••• 

r-oo foo roo 



poo poo poo _ 

+ dsi / dS2--- dS2n+lF{x + S2n+l)F{s2n+l + S2n)F{s2n + S2n~l) 

Jx Jx Jx 

■■■F{s2 + si)F{si + x) 
+ ■••, (B.7) 

_ _ 
Ki{x,x) = - dsiF{x + si)F{si+ x) 

J X 

roo roo roo _ _ 

+ dsi ds2 ds3F{x + S3)F{s3 + S2)F{s2 + si)F{si + x) 

J X J X Jx 



I X J X Jx 

+ ■■■ 

roo roo 



poo poo poo _ _ 

+ dSi / dS2--- dS2n+lF{x + S2n+l)F{s2n+l + S2n)F{s2n + S2n-l) 

J X J X J X 

■■■F{s2 + si)F{si + x) 
+ ■■■. (B-8) 

Appendix C: Proof of Proposition 4.1 

We remember F = ^2^^^ defined in §2^. The relation ( ^.231) with eq. ( 2.24 ) becomes 



= 2icr + - rQ("^)r, (c-i) 

for the cmKdV equations. We assume 



lim T = (ImC > 0), (C-2) 



and expand F as 



oo 



r-Ep^G,, (C.3) 



instead of eq. (|2.27|) . Substituting eq. ( C-3 ) for eq. we obtain a recursion formula, 



Gi+i = -5;,oi?^") + (Gi), + J2 G,Q^"''^Gi-„ / = 0, 1, • • • , (C-4) 
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where Q^"^^ and R^"^^ are given by 

2m-l 2m-l 

qM =vo1+ ^fc^fe' ^^"^ = -^oI+ E ^kek- (C-5) 

fe=l k=l 

We first show the following theorem. 

Theorem C. 1 Let X, Y, Z be 2"^"-^ x 2™"-"^ matrices given by 

2m- 1 

X = xol+ XkCk, (C-6a) 

k=l 

2 m- 1 

F = yoI+ E VkCk, (C-6b) 
fc=i 

2m- 1 

Z = 2:oI+ Y ^^^k' (C-6c) 
fc=l 

where the coefficients xq, a^^; yo) Vk'-, ^o-, Zk are real. Then there exist real numbers wq, Wk that 
satisfy 

2m-l 

W = XYZ + ZYX = wol+ Y ^kek- (C-7) 

k=l 

Proof By the use of eq. ( f4.2D , a direct calculation gives 

(2m-l \ / 2m-l \ / 2m-l 

i=l J \ j=l J \ k=l 

(2m-l \ / 2m-l \ / 2m-l 

k=l ) \ j=l ) \ i=l 

2m- 1 2m- 1 2m- 1 

= 2xoyo^oI + 2xoyo Y ^^^k + 2x0^0 ^ yjej + 2yozo ^ XiSi 

k=l j=l i=l 

{2m-l 1 f 2m-l f 2m-l 

X! Vj^kiejCk + ekCj) > + yos X! XiZkiciCk + 6^6^) > + 2;o< XiVji^^i^j + ej-ej 

j,fc=l J \i,k=l J [i,j=l 
2m- 1 

+ ^ ] XiyjZkis-iGje-k ~\~ S-k(ij^i) 
i,j,k=l 

2m- 1 2m- 1 2m- 1 

= 2xoyo^oI + 2xo?/o X! ^^^^^ + '^^ozo Y yj^i + ^^o^o X! ^^^^ 

fc=l j=l i=l 

2m-l 2m-l 2m-l 

-2x0 Y yiZi'^-2yo Y XiZil-2zQ ^ Xiyit 
1=1 1=1 1=1 

2m-l 2m-l 2m-l 2m-l 2m-l 2m-l 

i=l k=l k=l j=l j=l i=l 

This result shows that W does not include terms like CiCj or eiCjCk and can be expressed in the 
form of eq. ( |C^ ). □ 
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Using Theorem C. 1, we can prove by the inductive method that {Gi} are expressed by 

2m-l 
k=l 

where gj^^^ and ^f^^ are real coefficients. Therefore, T is given by 



2m~l 



oo (0) 2m-l / oo (A:) 

2m-l 

fc=i 

where ^^^\C) T^'^HC) satisfy 

We recall the asymptotic behavior of the Jost function at x — > ±00, 



/ 

O 



-i(x 



S(C)e'^^ 



as X ^ —00, 



as X —)■ +00. 



These relations yield 



lim ^le''^'' = I, lim ^W'^^ = A{C), 

x—>—oa - ■ ' -- 

where we have defined <pi and 02 by 



We easily see that 



lim (^atZ-r^ = (ImC > 0), 

(^1^00 



then we can replace T in eq. (p-10 ) with (f>2(f>i ^ ■ Thus det^(C) is expressed as 
det^(C) = exp{trlogyl(C)} 



expjtr J (log</>ie'^^)a;dx 
expJtr p {(t>ie'<^U(^ie'<^)-^dx 



(C-9) 



(C-10) 



(C-11) 

(C-12) 
(C-13) 

(C-14) 
(C-15) 
(C-16) 
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exp 



tr / Q('")</.2<Ar'dx 

^— oo 

exp|2™-i r (W°^(C)- E'W^(C))dxJ, (C-17) 



— oo 



fc=l 



where we have used eqs. ( [4.4[ ) and (4^). Due to eq. ( |C-11D , det^((^), as a function of C, satisfies 



det^(C) = {det^(-C*)}*. (C-18) 
This is the proof of Proposition 4.1 (1). Further, we obtain 
i?(C)A(C)-i= hin 020r'e-'^«" 

x— >+oo 



hm 

x— >+oo 



2m-l 



fc=i 

2m- 1 



^W(C)I+ E ^^'\C)ek, (C-19) 
fe=i 



with conditions 



r(0)(C) = {r(0)(-C*)}*, rW(C) = {rW(-C*)r- (C-20) 
Using eqs. (|C-ig| ) and (|C-20D , it is straightforward to prove Proposition 4.1 (2) (3). 
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Note added in proof — Equation ( 4.27 ) includes not only pure soliton solutions but also breather 
solutions. We should impose appropriate conditions on the arbitrary parameters in the residue 
matrices, e.g., C2j-iC2j = C2jC2j-i = C2jC2j-i = C2j-iC2j = O, or equivalently 

2m-l 



j=0 



to obtain pure soliton solutions. In this case, the 1-soliton solution for the cmKdV equations ( 2.55 ) 
is given by 

1 

2m-l \ 2 

Q^'-) (x, t) = 2r?i sech{2r?ix - 8r/i (r/f - 3^?)* - xq} ( 2 J] {cf^^^ 

\ 1=0 

where xo is defined by 

(2m-l \ ^2 

2 y: 
1=0 



For instance, vo{x,t) (cf. (|2.55| )) is given by 



2m-l 

(0|2 



voix,t) = 2r]i sech{277ix - 8r]i{r]j - 3£,l)t - xq} 2 ^ |c{ 

\ 1=0 

.{c(0)*e-2i?ix-8iCi(e2-3r?2)t ^ g(0)g2iCix+8iei(e?-3»?2)t| 
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